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Abstract. This paper is concerned with the existence of globally smooth solutions for 
the second boundary value problem for Monge- Ampere type equations and the application 
to regularity of potentials in optimal transportation. The cost functions satisfy a weak form 
of the condition A3, which was introduced in a recent paper with Xi-nan Ma in conjunction 
\Q | with interior regularity. Consequently they include the quadratic cost function case of 

Caffarelli and Urbas as well as the various examples in the earlier work. The approach is 
through the derivation of global estimates for second derivatives of solutions. 
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value problem for equations of Monge- Ampere type and its application to the regularity 
of potentials in optimal transportation problems with non-quadratic cost functions. 

X 
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This paper is concerned with the global regularity of solutions of the second boundary 



The Monge- Ampere equations under consideration have the general form 

(1.1) det{D 2 u - A{-, u, Du)} = B(-, u, Du), 

where A and B are given n x n matrix and scalar valued function defined onOxRx R" , 
where O is a domain in Euclidean n-space, W 1 . We use (x,z,p) to denote points in 
O x 1 x 1™ so that A(x,z,p) 6l"x R n , B(x,z,p) G R and (x,z,p) efixRx R n . The 
equation (1.1) will be elliptic, (degenerate elliptic), with respect to a solution u G C 2 {Vt) 
whenever 

(1.2) D 2 u- A(-,u,Du) > (>0), 
whence also B > (> 0). 
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A particular form of (1.1) arises from the prescription of the Jacobian determinant of 
a mapping Tu defined by 

(1.3) Tu = Y(-,u,Du), 

where Y is a given vector valued function on 11 x K x IR n , namely 

(1.4) det DY(-,u,Du) = ip(-,u,Du). 
Assuming that the matrix 

(1.5) Y p =[D Pi Y*] 

is non-singular, we may write (1.4) in the form (1.1), that is, 



(1.6) det{D 2 u + Y- 1 (Y x + Y z ®Du)} = 



\det Y p \ 



for degenerate elliptic solutions u. 

The second boundary value problem for equation (1.4) is to prescribe the image 

(1.7) Tu(Q) = O*, 

where Q* is a given domain in M n . When Y and i\) are independent of z and tfj is separable 
in the sense that 

(1.8) ip(x,p) = f(x)/goY(x,p) 

for positive /, g e L l {VL), L x (0*) respectively, then a necessary condition for the existence 
of an elliptic solution, for which the mapping T is a diffeomorphism, to the second 
boundary value problem (1.4) (1.7) is the mass balance condition 



(1-9) f f = f g. 



The second boundary value problem (1.4) (1.7) arises naturally in optimal transportation. 
Here we are given a cost function c : W 1 x W 1 — > R and the vector field Y is generated 
by the equation 

(1.10) c x (x,Y(x,p))=p, 

which we assume to be uniquely solvable for p e M n , with no n- vanishing determinant, 
that is 

(1.11) detc x , y (x,y) ^ 
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for all x,y G O x fi*. Using the notation 

d d 3d 



(1-12) Cij... jkl ... = 



dxi dxj dy k dyi 
we have 

(1.13) Y p (x,p) = [c i > j (x,Y(x,p))}, 

where [c l ' J ] is the inverse of [chj]. The corresponding Monge-Ampere equation can now 
be written as 

(1.14) det{D 2 u-c xx {-,Y(-,Du))} = |det c K)J/ |V>, 
that is in the form (1.1) with 

(1.15) A(x,z,p) = c xx (x,Y(x,p)), 

B(x,z,p) = \detc Xi y(x,Y(x,p))\ip(x,z,p). 

In the case of the (quadratic) cost function 

(1.16) c(x,y) = x-y, 
we have 

(1.17) Y(x,p)=p, Tu = Du, 

and equation (1.14) reduces to the standard Monge-Ampere equation 

(1.18) det D 2 u = i>. 

For this case global regularity of solutions was proved by Delanoe [D] , Caffarelli [C2] and 
Urbas [Ul], with (conditional) interior regularity shown earlier by Caffarelli [CI]. In this 
paper we will prove global estimates and regularity under corresponding conditions. In 
particular, we will assume that the cost function c G C 4 (M n x M n ) satisfies the following 
conditions: 

(Al) For each p, q G M n , there exists unique y = Y(x,p), x = X(q, y) such that 

c x (x,y)=p VxGO, 

c y (x,y) = q Vj/efi*. 



(A2) det c Xj y(x,y) ^ 0, V x G O, y G O . 
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(A3w) 



F(x,P;£,v) : = D PiP .A kl (x,p)£iZ j Ti k Tn 

>o V iGO^er^iieK". 



Conditions Al and A2 are the same conditions as in [MTW] but condition A3w is the 
degenerate form of condition A3 in [MTW] . As will be seen in our examples, we do not 
necessarily require c to be defined on all of lR n x R n and the vectors p and q in conditions 
Al and A3w need only lie in the ranges of c x (x, y) and c y (x, y) on O x O*. Moreover, as 
done at the outset in [MTW] , we may also write 

(1.19) F{x,p;£,ri) = (c ij>s - c k ' > 1 ' c ijik ,ci, :rs )c r ' k c s > 1 (x, y)^jr]kVi 

where y and p are related through Al. This shows that condition A3w is also symmetric 
in x and y. 

In our paper [MTW] , we also introduced a notion of convexity of domains with respect 
to cost functions, namely O is c-convex, with respect to O*, if the image c y (-,y)(Q) is 
convex in IR n for each y £ O*, while analogously O* is c*-convex, with respect to O, if the 
image c x (x, •)(&*) is convex for each x in O. For global regularity we need to strengthen 
these conditions in the same way that convexity is strengthened to uniform convexity. 
Namely we define O to be uniformly c-convex, with respect to O*, if O is c-convex, with 
respect to O*, <90 £ C 2 and there exists a positive constant 5 such that 

(1.20) [Dilj{x) - c i ' fc c lJ , ; (a:,y)7 fc (x)]r l r J (x) > 5 

for all x £ <90, y £ O*, unit tangent vector r and outer unit normal 7. By pulling back 
with the mappings c y (-,y), we see that this is equivalent to the condition that the image 
domains c y (-,y)(Q) be uniformly convex with respect to y £ O*. Similarly we call O* 
uniformly c*-convex, with respect to O, when c*(x, y) = c(y, x). Note that if O is simply 
connected with boundary 90 £ C 2 , then O is c— convex if and only if (1.20) holds for 
5 = 0. 

We can now formulate our main estimate. 

Theorem 1.1. Let c be a cost function satisfying hypotheses Al, A2, A3w, with respect 
to bounded C 4 domains 0,0* £ IR n which are respectively uniformly c-convex, c* -convex 
with respect to each other. Let ip be a positive function in C 2 (0 x 1R x M n ). Then any 
elliptic solution u £ C 3 (0) of the second boundary value problem (1.14), (1-7) satisfies 
the a priori estimate 

(1.21) \D 2 u\<C, 

where C depends on c,ip,Cl , O* and sup^ \u\. 
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As we will indicate later, the smoothness assumption on the solution and the data may 
be reduced. Further regularity also follows from the theory of linear elliptic equations 
for example if c, fi, are C°° then the solution u G C°°(0).The dependence of the 

estimate (1.21) on sup n \u\ may be removed if ip is independent of u as in (1.8). 

As a consequence of Theorem 1.1, we may conclude existence theorems for classical 
solutions by the method of continuity. 

Theorem 1.2. Suppose in addition to the above hypotheses that the function ip satisfies 
(1.8) (1.9). Then there exists a unique (up to additive constants) elliptic solution u G 
C 3 (Q) of the second boundary value problem (1.14), (1-V- 

From Theorem 1.2, we also obtain an existence result for classical solutions of the 
Monge-Kantorovich problem in optimal transportation. As above we let c G C 4 (IR n x R n ) 
be a cost function and Q, O* be two bounded domains in M. n satisfying the hypotheses 
of Theorem 1.1. Let / > 0, G C 2 (fi), g > 0, G C 2 (fi*) be positive densities satisfying the 
mass balance condition (1.9). Then the corresponding optimal transportation problem is 
to find a measure preserving mapping Tq : 0^0* which maximizes the cost functional 

(1.22) C(T)= [ f(x)c(x,T(x))dx 

Jn 

among all measure preserving mappings T from O to O*. A mapping T : O — > O* is 
called measure preserving if it is Borel measurable and for any Borel set £c!l*, 

(1-23) / f = [ g. 

Jt- 1 (e) Je 

Theorem 1.3. Under the above hypotheses, there exists a unique diffeomorphism T G 
[C 2 (0)] n maximizing the functional (1.22), given by 

(1.24) T(x) =Y(x,Du(x)), 

where u is an elliptic solution of the boundary problem (1.7), (1.14). 

The solution u of (1.7), (1.14) is called a potential. Note that in [MTW] and else- 
where the cost functions and potentials are the negatives of those here and the optimal 
transportation problem is written, (in its usual form), as a minimization problem. 

The plan of this paper is as follows. In Section 2, we prove that boundary conditions 
of the form (1.7) are oblique with respect to functions for which the Jacobian DT is non- 
singular and we estimate the obliqueness for solutions of the boundary value problem 
(1.14), (1.7) under hypotheses Al and A2, (Theorem 2.1). Here the twin assumptions 
of and fT being uniformly c and c*-convex with respect to each other are critical. 
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In Section 3, we prove that second derivatives of solutions of equation (1.14) can be 
estimated in terms of their boundary values under hypothesis A3w, (Theorem 3.1). This 
estimation is already immediate from [MTW] when the non-degenerate condition A3 
is satisfied. The argument is carried out for equations of the general form (1.1) (with 
symmetric A), in the presence of a global barrier, which is not necessary in the optimal 
transportation case (Theorem 3.2). This estimation also arises in the treatment of the 
classical Dirichlet problem [T3]. The proof of the global second derivative estimates 
in Theorem 1.1 for solutions of the boundary value problem (1.14), (1.7) is completed 
in Section 4. Here the procedure is similar to that in [LTU] and [Ulj.We remark here 
that this global estimate also extends to the more general prescribed Jacobian equation 
(1.6) [T3]. In Section 5, we commence the proof of the existence result, Theorem 1.2, 
by adapting the method of continuity [GT] and establish the result under a stronger 
uniform c-convexity hypothesis, (employed in earlier versions of this paper). Section 6 is 
devoted to the applications to optimal transportation and the proof of Theorem 1.3 from 
Theorem 1.2, which implies the global regularity of the potential functions in [MTW], 
under condition A3w. In Section 7, we finally complete the proof of Theorem 1.2 in its 
full generality, by showing that there exists a smooth function satisfying the ellipticity 
condition (1.2), together with the boundary condition (1.7), (at least for approximating 
domains). In the last section, we discuss our results in the light of examples, most of 
which are already given in [MTW]. 



2. Obliqueness 

In this section, we prove that the boundary condition (1.7) implies an oblique boundary 
condition and estimate the obliqueness. First we recall that a boundary condition of the 
form 

(2.1) G(x,u,Du) = on 90 

for a second order partial differential equation in a domain O is called oblique if 

(2.2) G P -7>0 

for all (x,z,p) G 90 x R x M n , where 7 denotes the unit outer normal to 90. Let us 
now assume that ip and 92* are C 2 defining functions for O and O* respectively, with 
V?, v?* < near 90, 90*, ip = on 90, tp* = on 90*, Vp, Vv?* ^ near 90, 90*. Then 
if u G C 2 (0) is an elliptic solution of the second boundary value problem (1.4), (1.7), we 
must have 



(2.2) 



cp*(Tu) = on 90, <p*(Tu)<0 near 90. 
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By tangential differentiation, we obtain 

(2.3) ^(DjT^t^O 

for all unit tangent vectors r, whence 

(2-4) ( p*(D j T i )= Xlj 

for some x > 0. Consequently 

(2.5) tfc l ' k (u jk -c jk ) =xiv 

that is 

(2-6) V*c l ' k w jk = xij, 

where 

(2.7) ^ij — ^ij ^ij '• 

At this point we observe that x > on <90 since |V<^*| ^ on <90 and det DT ^ 0. 
Using the ellipticity of (1.4) and letting [w^] denote the inverse matrix of [?%], we then 
have 

(2-8) tp*c i > k = X w jk lj . 

Now writing 

(2.9) G(x,p) = <p*(Y(x,p)), 
we have 

(2.10) p k :=G Pk (;Du) = X w jk lj , 
whence 

(2.11) /? • 7 = X«> ij 7i7j 

> 

on 90. We obtain a further formula for (3 ■ 7, from (2.6), namely 
(2-12) ^ k w jW tc l ' j = xVi&'Tj 
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Eliminating \ from (2.11) and (2.12), we have 

(2.13) ((3 • 7 ) 2 = (w'^aj^WMC^ 1 ^*). 

We call (2.13) a formula of Urbas type, as it was proved by Urbas [Ul] for the special 
case, c(x, y) = x • y, Y(^p) = p, of the Monge- Ampere equation. Note that to prove 

(2.13) , we only used conditions Al and A2 and moreover (2.13) continues to hold in the 
generality of (1.4). 

Our main task now is to estimate (3 ■ 7 from below for solutions of (1.14), (1.7). For 
this in addition to conditions Al, A2, we also need the uniform c and c* convexity of O 
and O* respectively. Our approach is similar to [Ul] for the special case of the Monge- 
Ampere equation and begins by invoking the key idea from [T] for estimating double 
normal derivatives of solution of the Dirichlet problem. Namely we fix a point xq on 
<90 where j3 • 7 is minimized, for an elliptic solution u G C 3 (0), and use a comparison 
argument to estimate 7 • D(f3 • 7) from above. Without some concavity condition in p 
the quantity (3 ■ 7 does not satisfy a nice differential inequality so we will get around this 
by considering instead the function 

(2.14) v = [3 -<y - nip*(Tu) 

for sufficiently large k, where now the defining function cp* is chosen so that 

(2.15) (D ij<P *(Tu) -c^c^ji-^Dk^iTu)^ > 5*|£| 2 

near 90, for all £ G M n and some positive constant 8q. Inequality (2.15) is possible by 
virtue of the uniform c*-convexity of O*, with the function if* given, for example by 

(2.16) v?* = ~ad* + b(d*) 2 , 

where a and b are positive constants and d* denotes the distance function for O*, [GT]. 
By differentiation of equation (1.14), in the form (1.1), we obtain, for r = 1, • • • ,n, 

(2.17) w ij {D ijUr - D Pk Aij(x,Du)D k u r - D Xr A ZJ (x, Du)} = D r \ogB. 
Introducing the linearized operator L, 

(2.18) Lv = w lJ (D ljV - D Pk AijD k v), 
we need to compute Lv for v given by (2.14). Setting 



(2.19) 



F(x,p) = G p (x,p) ■ 7(0:) - kG(x,p), 
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where G is defined by (2.9), we see that 



(2.20) v(x) = F(x,Du(x)). 
Writing b\ J = —D Pk Aij, we then have 

Lv =w lJ {F Pr DijU r + F PrPa D ir uD js u 

(2.21) + F XiX . + 2F XiPr D jr u + b k J (F Xk + F Pr D k u r )} 

In the ensuing calculations, we will often employ the following formulae, 

c l l l J (x,y)=D Xk c^(x,y) 

(2.22) =-c i ' l c r ^c kl , r (x,y), 

c; k J (x,y)=D yk c^(x,y) 

as well as (1.13). Indeed, using (1.13) and (2.22), we have 



(2-23) G PiPj =D Pj (rtc k n 



so that 



|2 



(2.24) G PtPj (x,Du)Uj >S* £ |c<%, 

for a further positive constant Kq. By choosing k sufficiently large, we can then ensure 
that 

(2.25) F PiP .(x,Du)te j <-±K\£\ 2 

near <90. Substituting into (2.20) and using (2.16), it follows that 

(2.26) Lv < --awn + C(w u + 1) + D Pk logBD k v, 

where C is a constant depending on c, ip, O and O*, as well as k. 

Next we observe that unless the defining function (p* extends to all of O* so that (2.15) 
is satisfied for all Tu 6 O*, we have no control on the neighbourhood of <90, where (2.26) 

9 



holds. This is remedied by replacing G in (2.19) by a function satisfying (2.24) in all of 
ft , agreeing with (2.9) near <9ft, for example by taking 



(2.27) G{x, P ) = m h {{^{Y{x^)), ai {\ P \ 2 -K 2 )}, 

where a\ and K are positive constants, with a\ sufficiently small and K >max|Z>u|, and 
for h sufficiently small, rrih is the mollification of the max-function of two variables. 

A suitable barrier is now provided by the uniform c-convexity of ft which implies, 
analogously to the case of ft* above, that there exists a defining function <p for ft satisfying 

(2.28) [Dij<p - ^CijjfaTvtDkdtej > 5 |£| 2 , 

in a fixed neighbourhood of <9ft, (for a constant do > 0). By appropriate choice, of say 
the constants a and b in (2.16), without the *, (or following the uniformly convex case 
in [GT, Chapter 14]), we may obtain, by virtue of (2.21), 

(2.29) Lip > 5 w li + Kw ij DiipDjip, 

for a given constant K. Combining (2.26) and (2.29), and using the positivity of B, we 
then infer by the usual barrier argument, ( which entails fixing a small enough neighbour- 
hood of 0ft, [GT]), 

(2.30) -y-Dv(x )<C, 

where again C is a constant depending on c, ft, ft* and ip. From (2.28) and since xq is a 
minimum point of v on 9ft, we can write 

(2.31) Dv(x ) = T-f(x ) 

where r < C. To use the information embodied in (2.31), we need to calculate 
(2.32) 

=V^A(T^)c^ 7j + ipl^Da, + <PllM' j + cfD^u) 

=ip k C '■'{Dilj - C ' Cij^r) + {ip kl - ip r C ' C S)k l)C ,J ^jDiT u 

Multiplying by (/?*c*' 1 and summing over i, we obtain 

(2.33) ri^'DiiP ■ 7 ) =vfe*c fc 'V'*(A7, - c s > r Cij , slr ) 

+ tfc^C^^jC^Wimiipli - ip*C r ' S C s ,kl) 

>5o5>!c^'| 2 
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by virtue of the uniform c-convexity of O, the c*-convexity of O* and (2.6). Consequently, 
from (2.19) and (2.31), we obtain at xq, 

(2.34) -KWu&WipWj < C((3 • 7) - r 

for positive constants, C and tq. Hence if j3 • 7 < ro/2C, we have the lower bound 

(2.35) wuWtfVj > 



To complete the estimation of f3 ■ 7 we invoke the dual problem to estimate w lJ ; 7*7,7 at 
xq. Assuming for the moment that Tu is one to one, we let u* denote the c-transform of 
u, defined for y = Tu{x) G fi* by 

(2.36) u*(y) = c(x, y) - u(x). 
It follows that 

(2.37) Du*(y) = cy(x,y) 

= c y (T*u*(y),y), 

where 

(2.38) T*u*(y) = X(Du*,y) 

= {Tu)-\y\ 

and the second boundary value problem (1.14), (1.7) is equivalent to 

(2.39) \detD y (T*u*)\=g(y)/f(T*u*) in Q*, 

(2.40) T*0* = n. 

Noting that the defining functions ip and ip* may be chosen so that V(p = 7, V</?* = 7* 
on 90, <9Q* respectively, we clearly have for x G 90, j/ G Tu{x) G 90*, 

(2.41) )9-7(aj) = c*' 4 (^yWfc(y) 

= ^-7*(y), 

where 



(2.42) 



^(y) = %(r(V,!/))- 
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Hence the quantity (3* ■ 7* is minimized on <90* at the point y = Tu(x ). Furthermore, 
for y = Tu{x), x G 90, 

(2.43) w ij 7aj(x) = wli{y)c Ki c l ^{x,y)ip iVj {x), 

where 

(2-44) w* kl (y) = u* ykyi (y)-c kl (x,y). 

Applying now the estimate (2.35) to u* at the point yo G 90*, we finally conclude from 

(2.13) the desired obliqueness estimate 

(2.45) 0-7 > 5 

on 90 for some positive constant 5 depending only on O, O*, c, and ip. 

The above argument clearly extends to arbitrary positive terms B (1.15). Noting also 
that it suffices in the above argument that T need only be one-to-one from a neighbour- 
hood of the point xq to a neighbourhood of yo, we have the following theorem. 

Theorem 2.1. Let c G C 3 (lR n xR n ) be a cost function satisfying hypotheses Al, A2, with 
respect to bounded C 3 domains 0,0* C M n , which are respectively uniformly c-convex, 
c* -convex with respect to each other. Let if) be a positive function in C 1 (0 x H. x R n ). 
Then any elliptic solution u G C 3 (0), of the second boundary value problem (1.14), (1-7) 
satisfies the obliqueness estimate (2.45). 

Note that Tu is automatically globally one-to-one under the hypotheses of Theorem 
1.2 by virtue of the change of variables formula. In ensuing work , (see [T3]), we extend 
Theorem 2.1 to the more general prescribed Jacobian equation (1.4). The main difference 
is that we cannot directly use the c-transform to get the complementary estimate to (2.35) 
Instead the quantities there are transformed using the local diffeomorphism Tu. Indeed 
we could also have avoided the use of duality in the proof of Theorem 2.1 by direct 
transformation of (2.31). 

3. Global second derivative bounds 

In this section we show that the second derivatives of elliptic solutions of equation 

(1.14) may be estimated in terms of their boundary values. For this estimation and the 
boundary estimates in the next section, it suffices to consider the general form (1.1) under 
the assumption that the matrix valued function A G C 2 (0 x R x M n ) satisfies condition 
A3w, that is 



(3.1) 
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for all (i,z,j))6llxRx IR n , £, rj G IR n , £ _L rj. We also assume A is symmetric, which is 
the case for the optimal transportation equation (1.14). When (3.1) is strengthened to 
the condition A3 in [MTW], that is 

(3-2) D PkPl A i:j (x,z,p)^jVkVi > <5|£| 2 M 2 

for some constant 6 > 0, for all (i,z,p) GOxRx IR n , £, 77 G R n , £ _L 77, then the global 
second derivative estimate follows immediately from our derivation of interior estimates 
in [MTW]. In the general case the proof is much more complicated and we need to also 
assume some kind of barrier condition, namely that there exists a function (p G C 2 (0) 
satisfying 

(3.3) [Dij<p(x) - D Pk A ij (x,z,p)D k (p(x)]^ j > 5|£| 2 

for some positive 5 > and for all £ G M n , x, z,p G some set f/ C O x 1 x M n , whose 
projection on Q is O. In general, condition (3.1) implies some restriction on the domain 
Q, but for the case of equations arising in optimal transportation, it can be avoided by 
a duality argument. 

Our reduction to the boundary estimation follows the approach in [GT], originating 
with Pogorelov, with some modification analogous to that in [LTU]. Let u G C 4 (0) be 
an elliptic solution of equation (1.1), with x, u(x), Du(x) G U for x G O and £ a unit 
vector in R n . Let v be the auxiliary function given by 

(3.4) v = v(-, £) = logiwij&Q + t\Du\ 2 + kv, 
where Wij = DijU — Aij. By differentiation of equation (1.1), we have 

(3.5) w lJ [DijU^ - D^A i:j - (D z Aij)u£ - {D Pk Aij)D k u^\ 

=D^B + (D Z B) U £ + (D Pk B)D kU £, 

where B = logS. A further differentiation yields 
(3.6) 

w tJ [DijUtf - D^Aij - (D zz A ij )(u^) 2 - {D PkPl Aij)D k u^Diu^ 

- (D z A i:j )u^ - (D Pk Aij)D k u^ - 2(D{. z A ij )u{ i 

- 2{D^ Pk A ij )D k u i - 2(D zpk Aij)(D k u^u 6 ] - w ik w jl D^ Wij D c w k i 
= D#B + {D zz B)u\ + (D PkPl B)D kH D lUc 

+ 2{D iz B) Ui + 2{D iPk B)D kUi + 2{D zpk B){D k u i )u i 

+ (D z B)u^ + (D Pk B)D k u^. 
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Furthermore differentiating (3.4) we have 



(3.7) 

DiV = + 2rD k uD ki u + kD^, 

(3-8) 

DijV = — — — 2 + 2T(D ik uD jk u + D k uD ijk u) + nD^y. 

where we have written = Dijw^j. Using condition A3w in (3.6) and retaining all 
terms involving third derivatives, we estimate 

(3.9) Lu^ : = w ij (DijUtf + b% D k u^) - (D Pk B)D k u^ 

> w^w^D^WijD^WH - C{{1 + w u )w u + (w u ) 2 } 

where, as in the previous section, b\ J = —D Pk Aij and C is a constant depending on the 
first and second derivatives of A and logS and sup n (|u| + \Du\). To apply A3w, we fix 
a point x G O and choose coordinate vectors as the eigenfunctions of the matrix [wij] 
corresponding to eigenvalues < Ai < • • • < X n . Writing Aij^i = D PkPl Aij, we then 
estimate 

w tJ A ijM u k ^u^ > w tJ A ijM w k ^w^ - Cw"(l + wu) 

> Y A rr,kl( X kCk)(XlCl)-Cw U (l + W U ) 
k or l=r T 

> -C{w ii (l + w ii ) + wu} 



From (3.9), we obtain also 

(3.10) Lwk > w tk w jl D^ WlJ D^w kl - C{(1 + w u )w u + w u } 



for a further constant C. Here we use equation (3.5) to control the third derivative 
term arising from differentiating A k i£ k £i. From (3.8) and (3.10), we obtain, after some 
reduction, 

1 1 

(3.11) Lv > w lk w jl DtwaD £ Wki ttW 13 DiWctDiWct 

;r ii air 

+ 2tw u + kw u - C{ [(1 + w u )w 11 + w u ] + t + k] 

Now suppose v takes its maximum at a point xq G O and a vector £, which we take to 

be e\. We need to control the first two terms on the right hand side of (3.11). To do this 
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we choose remaining coordinates so that [wij] is diagonal at xq. Then we estimate 

1 1 

(3.12) w lk w jl DfWiiDtw k i n-w 13 DiWftDiWee 

ICrr llir 

- 1 -w ii w^{D 1 w ij ) 2 - -L w "(D lWll ) 2 



Wn w 
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> _i_ ^^"(D^?. _ w ii {D i w 1 if] 
w ii i>i 

11 i>i 

2 ^ 

H 2~ 2^^ M [-Diwii - A^ii] [-Di^h + DiWu] 

W 1 i 

11 2>1 

i>i 

w n i>i 

> -Cw" 

Combining (3.11) with (3.12), we obtain the estimate, at xo, 

(3.13) Lv > tw u + kw" - C{t + k}, 

for either r or k sufficiently large. Note that when we use (3.7) in the second last line of 
(3.12), we improve (3.13) by retention of the term ^2 i>1 w ll (DiWn) 2 on the right hand 
side , which corresponds to the key term in the Pogorelov argument for interior estimates 
[GT]. 

From (3.13), we finally obtain an estimate from above for wa(xo), which we formulate 
in the following theorem. 

Theorem 3.1. Let u G C 4 (0) be an elliptic solution of equation (1.1) in , with 
x,u(x),Du(x) G U, for all ifO. Suppose the conditions A3w and (3.3) hold and B is 
a positive function in C 2 (0 x R x W 1 ). Then we have the estimate 

(3.14) sup|D 2 w| < C(l + sup|D 2 w|), 

n an 

where the constant C depends on A, B, and s\ipu(\z\ + \p\). 

Note that we only need the condition A3w to hold on the set U. 

From the proof of Theorem 3.1 we obtain the corresponding estimate for equation 
(1.14), without the barrier condition (3.3). 
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Theorem 3.2. Let u G C 4 (0) be an elliptic solution of equation (1.14) ^ n ^ with 
Tu{Vl) C O*. Suppose the cost function c satisfies hypotheses Al, A2, A3w and B is 
a positive function in C 2 (0 x R x R n ). Then we have the estimate (3.14). 

To prove Theorem 3.2, we take k = in the proof of Theorem 3.1, to obtain an 
estimate for wu in terms of w l \ that is 

(3.15) wu <esupw u + C e (l + sup \D 2 u\), 

for arbitrary e > 0, with constant C £ also depending on e. If T is globally one-to-one, 
we then conclude (3.13), in the optimal transportation case, by using the dual problem 
(2.37), (2.38). More generally, we consider the dual function v* in place of (3.4), given 
by 

(3.16) = v{x,£) = \og(w^ c^Cj^i) + r\c y (x,Tu(x))\ 2 

and suppose it is maximized at a point Xq in O. Since T will now be one-to-one from a 
neighborhood Af of Xq to a neighbourhood Af* of = Tu(xq), we may then proceed as 
before, noting that in Af*, v* is given by (3.4) with u replaced by its c— transform u* . 

The estimate (3.14) arose from our investigation of the classical Dirichlet problem , 
(see [T3] ). We remark also that from (3.16), we see that (3.3) is also not needed when 
n = 2. 

4. Boundary estimates for the second derivatives 

This part of our argument is similar to the treatment of the oblique boundary value 
problems for Monge-Ampere equations in [ITU, U3]. The paper [LTU] concerned the 
Neumann problem, utilizing a delicate argument which did not extend to other linear 
oblique boundary conditions. For nonlinear oblique conditions of the form (2.1) where 
the function G is uniformly convex in the gradient, the twice tangential differentiation of 
(2.1) yields quadratic terms in second derivatives which compensate for the deviation of 
[3 — G p from the geometric normal and permit some technical simplification for general 
inhomogeneous terms ib [U3]. 

First we deal with the non-tangential second derivatives. Letting F e C 2 (Q x R x R n ) 
and v = F(-,u, Du), where u G C 3 (0) is an elliptic solution of equation (1.1), we have 
from our calculation in Section 2, 

(4.1) \Lv\ <C(w" + w u + 1), 
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where L is given by (2.17) and C is a constant depending on A, B, G, O and |w|i ; n. Now- 
using the equation (1.1) itself, we may estimate 

(4.2) wf 1 < Cw«, 
so that, writing M = sup n wa, we have from (4.1) 

71-2 

(4.3) \Lv\ < C(l + M) — w 11 . 

Hence, if there exists a C 2 defining function (p satisfying (3.3) near 90, together with 
ip = on 90, we obtain by the usual barrier argument, taking F = G, 

(4.4) \D(0-Du)\<C(l + M)%=* 
on 90, so that in particular 

(4.5) W{3p<C(l + M) — 
on 90. Now for any vector £ G M n , we have 

(4.6) = w TT + b(w Tf 3 + W/3 T ) + b 2 wpp, 
where 

(4.7) b =7T L > T = t-b$. 

Suppose takes its maximum over 90 and tangential £, |£| = 1 at xq G 90 and £ = ei. 
Then from (4.5) and (4.6) and tangential differentiation of the boundary condition (2.1) 
we have on 90, 

(4.8) tun < |ei-6/3| 2 wii(0)+6F(-,w, J Dw) + C6 2 (l + M) — , 

for a given function F G C 2 (0 x R x R n ). Combining (2.26), (3.10), (4.1) and (4.2) 
and utilizing a similar barrier argument to that is Section 2, we thus obtain the third 
derivative estimate 

2ra-3 

(4.9) -Dp Wll (xo) <C(1 + M) — , 

Differentiating (2.1) twice in a tangential direction r, with t(xq) = ei, we obtain at xq, 

(4.10) (D PkPl G)u lk uu + (D Pk G)u llk < C(l + M), 
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whence we conclude from (4.9) 

o o 2 n — 3 

(4.11) max|L> 2 <u| < C(l +sup D 2 u ^ 

by virtue of the uniform convexity of G with respect to p. Taking account of the global 
estimate (3.13), we complete the proof of Theorem 1.1. □ 

Once the second derivatives are bounded, the equation (1.1) is effectively uniformly 
elliptic so that from the obliqueness estimate (2.43), we obtain global C 2,a: estimates 
from the theory of oblique boundary value problems for uniformly elliptic equations in 
[LT] . By the theory of linear elliptic equations with oblique boundary conditions [GT] , 
we then infer estimates in C 3 ' a (0) for any a < 1 from the assumed smoothness of our 
data. We may also have assumed that our solution u G C 2 (0). 

As in the previous section, the technicalities are simpler when condition A3w is 
strengthened to condition A3 and we also obtain local boundary estimates for the second 
derivatives. To see this we estimate the tangential second derivatives first by differenti- 
ating the equation (1.1) and boundary condition (2.1) twice with respect to a tangential 
vector field r near a point y G dfl. We then obtain an estimate for r]D TT u, for an appro- 
priately chosen cut-off function r\. The mixed tangential-normal second derivatives D Tn u 
are estimated as above by a single tangential differentiation of (2.1) so that the double 
normal derivative may be obtained either from (4.5) or from the equation (1.1) itself and 
the estimates in Section 2 for w %3 ^i r )j from below, similarly to the Dirichlet problem, see 
[Tl]. 

5. Method of continuity 

To complete the proof of Theorem 1.2, we adapt the method of continuity for nonlinear 
oblique boundary value problems, presented in [GT] and already used in the special case 
(1.16) (1.17) [Ul]. The situation here is more complicated unless we know in advance 
that there exists a smooth function satisfying the ellipticity condition (1.2) together 
with the boundary condition (1.7). Later in Section 7, we shall prove the existence of such 
a function, (at least for approximating domains). Otherwise we need to consider families 
of subdomains. To commence the procedure, we fix a point xq G O. Then for sufficiently 
small radius r > 0, the ball fi = B r (xo) C O will be uniformly c-convex with respect to 
fi* and the function uq, given by 

(5.1) uq{x) = - \x - x \ 2 +p ■ (x - x ), 

will satisfy the ellipticity condition (1.2). Moreover the image Oq = Tu (Q ) will be 
uniformly c*-convex with respect to O with Tu a diffeomorphism from fi to Oq- To 

18 



see this we observe that 



(5.2) c x (x ,Oo) = B Kr (p Q ), 

so that by taking kt small enough, we can fulfill condition (1.19) on <90q, with respect 
to xo G O, for constant 5o = as large as we wish. Suppose now we can foliate O — O 
and O* — fig by boundaries of c-convex and c*- convex domains, respectively. That is 
there exist increasing families of domains {O t }, {O*}, < t < 1, continuously depending 
on the parameter t, such that 

(i) n t c o, o* c o*, 

(ii) Oi = 0,0* = O*, 

(iii) <90 t , <90^ G C 4 , uniformly with respect to t, 

(iv) 0^,0* are uniformly c-convex, c*-convex with respect to 0*,0, respectively. 
The construction of such a family is discussed below. 

Given our families of domains O t ,0^, < t < 1, we need to define corresponding 
equations. Let B be a positive function in C 2 (0 x W 1 ) and / a positive function in 
C 2 (0) such that 

(5.3) / = -au + log[det{D 2 w - c xx (-, Y(-, Du )}/B(-, Du Q )\ 

in O , for some fixed constant a > 0. We then consider the family of boundary value 
problems: 

(5.4) F[u] : = det{D 2 u-c xx (-,Y(-,Du)} = e au+{1 ~ t)s B(-,Du), 

Tu(n t ) = Y(;Du)(n t ) = n* t , 

From our construction and the obliqueness, we see that uq is the unique elliptic solution 
of (5.4) at t = 0. 

From Section 2, we also see that the boundary condition in (5.4) is equivalent to the 
oblique condition 

(5.5) G t {-, Du) := (pt(Y(-, Du)) = ondO t . 

To adapt the method of continuity from [GT], we fix a G (0, 1) and let £ denote the subset 
of [0, 1] for which the problem (5.4) is solvable for an elliptic solution u = u t G C 2 ' a (O t ), 
with Tu invertible. We then need to show that E is both closed and open in [0, 1]. First 
we note that the boundary condition (5.4) implies a uniform bound for Dut. Integrating 
the equation (5.4), we then obtain uniform bounds for the quantities 



e aUt . 
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so that the solutions u t will be uniformly bounded for a > 0. Uniform estimates in 
C 2,1 (Cl) now follow from our a priori estimates in Section 4, which are also clearly in- 
dependent of t G [0, 1]. By compactness, we then infer that £ is closed. To show E is 
open, we use the implicit function theorem and the linear theory of oblique boundary 
value problems, as in [GT]. The varying domains {fit} may be handled by means of 
diffeomorphisms approximating the identity, which transfer the problem (5.4) for t close 
to some to ^ E to a problem in Clt - We then conclude the solvability of (5.4) for all 
t G [0, 1], which implies there exists a unique elliptic solution u = u a G C 3 (0) of the 
boundary value problem 

(5.6) F[u] = e au B(-,Du), 

Tu(Cl) = Cl* 

for arbitrary a > 0, with Tu one-to-one. To complete the proof of Theorem 1.2, (at least 
when the above foliations exist), we assume that B satisfies (1.5), (1.8) and (1.9). As 
above we see that the integrals 




are uniformly bounded, with D(au a ) — > as a — > 0. Consequently au a — * constant = 
by (1.19) and modulo additive constants, u a — > u as a — > 0, where u is the solution of 
(1.14), (1.7), as required. 

We may construct the family of domains {Ot} used above, if we are given a defining 
function C defining function ip, satisfying 



(5.7) [Dij<p(x) - c l ' k Cij ,i(x, y)D k (p(x)]te j > 5 |^| 2 

for all x G O, y G CI*, £ G M. n , which takes its minimum at Xq. Note that the uniform 
c-convexity of CI implies the existence of a defining function satisfying (5.7) in a neigh- 
bourhood of dCl, as in (2.26). There are various ways of constructing suitable families 
from a global defining function, <p. In particular taking (p(xo) = —1, we may choose 

Cl t = {xeCl\ (p t {x) < 0} 

where <p t is defined by 



(5.8) (ft = t((p - a) + (t - t)cpo, 

for (po(x) = \x — x \ 2 — r 2 , t < t , for some < t < 1 and a close to —1 to first 

deform to a small sub-level set of (p, followed by taking <p t = (1 — t)a/ (1 — to) for t > t . 
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Alternatively, we could have chosen fi = ip < a at the outset and only used the second 
deformation.The domains Q£ may be similarly constructed. If the curvatures of <9fi* 
are sufficiently large, for example if fi* is a small ball, then the existence of a defining 
function satisfying (5.7) follows by pulling back from a single image c y (-, yo)(Q). As a 
byproduct, we see that if O is uniformly c-convex with respect to a single point yo G O*, 
there at least exists a smooth function uq G C 3 (0) satisfying the ellipticity condition 
(1.2). Moreover Tuq is a diffeomorphism from O to a small ball B r (yo) C O*. 

From the above considerations, we see that the proof of Theorem 1.2 is completed 
in the cases where either Q or Q* is a small ball. The general case will then follow by 
further use of the method of continuity if there exists a defining function satisfying (5.7) 
for either domain. However we will take up a different approach from this point and 
use the function uq constructed above to construct a further function u\ approximately 
satisfying our given boundary conditions to which the method of continuity can be applied 
without domain variation. Specifically we will prove 

Lemma 5.1. Let the domains O and fi* and cost function c satisfy the hypotheses of 
Theorem 1.1. Then for any e > 0, there exists a uniformly c* -convex approximating 
domain, O*, lying within distance e of Q* , and satisfying the corresponding condition 
(1.19) for fixed So, together with a function u\ G C 4 (0) satisfying the ellipticity condition 
(1.2) and the boundary condition (1.7) for Q* . 

From Lemma 5.1 we complete the proof of Theorem 1.2. We defer the proof of Lemma 
5.1 to Section 7 as the proof will use some of the same ingredients as in our discussion 
of optimal transportation in the next section. In Section 7, we will also indicate an 
alternative and direct construction of the function u±, which avoids domain variation 
altogether in the method of continuity. 



6. Optimal Transportation 

The interior regularity of solutions to the optimal transportation problem is considered 

in [MTW], under conditions Al, A2, A3 and the c*- convexity of the target domain O*. 

Our approach is to first show that the Kantorovich potentials are generalized solutions 

of the boundary value problem (1.14), (1.7) in the sense of Aleksandrov and Bakel'man. 

The c*-convexity of O* is used to show the image of the generalized normal mapping 

lies in O* and condition A3 is employed to obtain a priori second derivative estimates 

from which the desired regularity follows. The potential functions u and v solve the dual 
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problem of minimizing the functional 



(6.1) 



I(u,v) = fu + 



gv 




K = {(u,v) | u,v G C°(fi),C°(fi*) resp. u(x)+v(y) > c(x,y) for all x G O, y G fi*}. 
The potential functions (w, v) satisfy the relations 



that is they are the c* and c transforms of each other. Since c G C 1 ' 1 , they will be 
semi-convex. The optimal mapping T is then given almost everywhere by (1.23) and the 
equation (1.14) will be satisfied with elliptic solution u almost everywhere in O. The 
functions u and v are respectively c and c*- convex. A function u G C°(Q) is called 
c — convex in O if for each xq G O, there exists G M n such that 

(6.4) u(x) > u(x Q ) + c(x, y ) - c(x , y Q ) 

for all x G O. If w is a c— convex function, for which the mapping T given by (1.23) 
is measure preserving, then it follows that u is a potential and again T is the unique 
optimal mapping. These results hold under the hypotheses Al and A2 and it suffices to 
assume the densities /, g > 0, G L 1 (fi), L 1 (Q*), respectively, whence the mapping T is 
only determined almost everywhere on the set where f is positive. The reader is referred 
to [C3,GM,MTW,U2,V] for further details. 

From the above discussion we see that the solution of the boundary value problem 
(1.14), (1.7) will automatically furnish a potential for the optimal transportation problem 
if it is c-convex. Note that ellipticity only implies that the solution is locally c-convex and 
we need a further argument to conclude the global property, unlike the case of quadratic 
cost functions and convex solutions. First we recall the concept of generalized solution 
introduced in [MTW]. Let u be a c-convex function on the domain O. The c — normal 
mapping, \ u , is defined by 

(6.5) Xu(x ) = {y G R n | u(x) > u(x ) + c(x, y) - c(x , y ), for all x G O}. 

Clearly, Xu(%o) G Y(x , du(x )) where d denotes the subgradient of u. For g > 0, G 
L^ oc (IR n ), the generalized Monge- Ampere measure n[u,g] is then defined by 



(6.3) 



u(x) = sup {c(x, y) — v 
yen 



■(*)}, 



v(y) = sup {c(x, y)-u 



(6.6) 



^Kf](e) = / 9 

"'X«(e) 

22 



for any Borel set e G O, so that u satisfies equation (1.14) in the generalized sense if 



(6.7) fJ>[u,g] = fdx. 

The boundary condition (1.7) is satisfied in the generalized sense if 

(6.8) 0*Cx„(0), \{x G O I f(x) > and X u(x) - O V 0}| = 

The theory of generalized solutions replicates that for the convex case, c(x, y) = x ■ y, 
[MTW]. If / and g are positive, bounded measurable functions on 0,0* respectively 
satisfying the mass balance condition (1.9), and c satisfies Al, A2, then there exists a 
unique (up to constants) generalized solution of (6.7), (6.8), (with g = outside O*), 
which together with its c transform v, given by (6.3), uniquely solves the dual problem 
(6.1), (6.2)([MTW]). 

Now let u G C 2 (0) be an elliptic solution of the boundary value problem (1.7), (1.14) 
and v a c— convex solution of the corresponding generalized problem. By adding con- 
stants, we may assume inf^u — v) = 0. We need to prove u = v in O, that is the strong 
comparison principle holds. Let O' denote the subset of O where u > v and first suppose 
that dVt' n O 7^ 0. Note that if v G C 2 (0), this situation is immediately ruled out by 
the classical strong maximum principle [GT]. Otherwise we may follow the proof of the 
strong maximum principle as there will exist a point xq G 90' n O, where O' satisfies an 
interior sphere condition, that is there exists a ball Bcfi-O' such that xq G 90' n dB, 
u(xo) = v(xo) and u > v in B. Since v is semi-convex, v will be twice differentiate at 
xq, with Dv(xo) = Du(xo). Moreover by passing to a smaller ball if necessary we may 
assume both u and v are c— convex in B. Since u is a smooth elliptic solution of (1.14), 
there will exist a strict supersolution w G C 2 (B — B p ), for some concentric ball B p of 
radius p < R, satisfying w{xq) = u(xo), w > v on dB U dB p , Dw(xo) ^ Du{xq). By 
the comparison principle, [MTW], Lemma 5.2, we have w > v in B — B p , and hence 
Dw(xq) = Dv(xq), which is a contradiction. Thus we may assume <90 ; n O = 0, that 
is u > v in O with u(xq) = v(xq) for some point xq G 90. From our argument above, 
we obtain a function w G C 2 (B — B p ) satisfying w(xo) = u(xo) = v(xq), v < w < u in 
B — B p , together with 

(6.9) u(x) — w(x) > e\x — xo\ 

for all x G Br — B p . Since v < w in Br — B p , this contradicts the obliqueness condition 
(2.43) if O* is c* -convex. 

Alternatively we may proceed directly ( and more simply) as follows to show that the 
solution u is c-convex, using the property that Tu is one-to-one. In fact, as mentioned 
previously in Section 2, this would follow automatically from the change of variables 
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formula by virtue of the mass balance condition (1.9). Let xq G Vt and yo = Tu(xq). 
Suppose there exists a point x\ G O, where 



(6.10) 



u{x\) < c(x 1 ,y ) - c(x , y ). 



By downwards vertical translation, there exists a point x 2 G <90, satisfying 



(6.11) 



u(x) > u(x 2 ) + c(x, y ) - c(x 2 , y ). 



for all iGO. Putting y 2 = Tu(x 2 ) : we must also have 



(6.12) 



c x (x 2 ,y 2 ).^(x 2 ) < c x (x 2 ,y ).'j(x 2 ), 



which again contradicts the c*— convexity of O*. 
Remarks. 

In the first proof above, we employed a comparison result that if u is a classical 
elliptic supersolution of (1.14) dominating a generalized subsolution v on the boundary 
of a subdomain Q' , then u > v in Q' . In our local uniqueness argument in [MTW], we also 
used implicitly the complementing result that if u is an elliptic subsolution dominated by 
a generalized supersolution v on <9Q', then u < v in Q' . However, in this case, we cannot 
apply Lemma 5.2 in [MTW] directly as local c-convexity of v may not imply global c- 
convexity in O, unless v G C 1 (0). This situation is rectified in [TW], under the strong 
A3 hypothesis. However if O and O* lie respectively in domains fio and Qq satisfying 
the hypotheses of Theoreml.2, with O* also c*-convex with respect to fi , and f and g 
are positive in L 1 (0) and L 1 (0*) respectively, then it follows directly by approximation 
from Theorem 1.3 that the local c-convexity of the potential u solving the Kantorovich 
dual problem implies its global c-convexity. Other results, such as the c-convexity of 
the contact set under condition A3w, also follow from Theorem 1.3 by approximation. 
The reader is referred to Loeper [L] for a full treatment of this approach, including the 
sharpness of condition A3w for regularity. 



In this section, we provide the proof of Lemma 5.1, thereby completing that of Theorem 
1.2. For this purpose, we need to draw on a geometric property of c-convex domains 
introduced in [TW]. Namely, suppose that O is uniformly c-convex, with respect to O*, 
and that the cost function c satisfies conditions Al, A2 and A3w. Denoting as before 
the unit outer normal to 90 by 7, we see that the level set £ of the function e, given by 



7. Completion of Proof of Theorem 1.2 



(7.1) 



e(x) = e y (x) = c(x, y) - c(x, y ), 
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passing through xq, is tangential to dQ at xq if 

(7.2) y = y(x ,p +t7o), 

for t > 0, po = c x (xo, yo), 70 = 7(^0) that is, y lies on the c*-segment which is the image 
under Y(xq,-) of the line from po with slope 70. Then it follows from [TW] that O lies 
strictly on one side of £, whence 

(7.3) c(x, y) - c(x ,y) < c(x, y ) - c(x ,y ), x G Cl. 

To prove (7.3) directly from (1.19), we take x = 0, set x' = (xi, x n -i) and choose 

coordinates so that 70 = (0, , —1). By Taylor's formula, 

e(x) - e(x ) < -tx n + [Aij(0,p + £70) - Aij(0,p o )]x'ix' j + tO(\x\\x n \ + \x\ 3 ). 

Using (1.13), (1.19), condition A3w and again Taylor's formula, we have 

[Aij(0,p o +^70) - Aij(0,po)]x'ix'j < -tc l ' n c ijt i(xo,yo)x' ix' j 

< tD i 'Yj(xo)x' i x'j - tS \x'\ 2 , 

so that, 

e(x) — e(xo) < —tx n + tD i 'jj(xo)x'ix' j + tO(|x||x n | + |x| 3 ) < 0, 

for 16O - xq, sufficiently small. Consequently O lies locally, strictly on one side of S. 
We can then verify the global inequality (7.3) by contradiction, as in [TW]. For if (7.3) 
is violated, the set 

U a = {x G 90; e(x) > a\} 

contains two disjoint components, for sufficiently small a > 0. Increasing a, we see that 
these components will meet first at a point x* G <90 at which the level set of the function 
e is tangential, contradicting the local inequality at x*. 

Now, to commence the proof of Lemma 5.1, we take uq to be a function as constructed 
in Section 6, that is uq is a smooth uniformly c-convex function on O, whose c-normal 
mapping Tuq = Y(-,Duq) has image u* , which is a c*-convex subdomain of O*. Here 
we call a c-convex function uniformly c-convex if it also satisfies the ellipticity condition 
(1.3). We remark that the c-convexity of uq could also have been proved from (7.3), 
using the c-convexity of Q and condition A3w, instead of the c*-convexity of O* which 
we used in Section 6. Also by approximation, we may assume uq G C°°(0). A function 
h is called a c-function if it has the form 



K) = c (-iVo) + a 
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for y G M n and some constant a. When c(x, y) = x • y, a c-function is a linear function. 
Obviously the c-normal mapping of h is the constant map, Th(x) = yo for all x G M n . 
Let 

(7.4) ui(x) = sup{it , h(x)}, xeQ s , 

where the sup is taken over the set S of c-functions h with h < uq in O, 7^(0) C O* and 
O 6 = {x G K n : dist(x, O) < 5}, for some 5 > 0, is a neighbourhood of O. The following 
lemma describes the properties of u±. 

Lemma 7.1. . Assume that the cost function c satisfies Al, A2, A3w and the domains 
O and O* are uniformly c-convex with respect to each other. Then, for sufficiently small 
5, the function u\ is a c-convex extension of uq to O 5 , whose c-normal image under u\ 
is O . Moreover for any point igO { -!1, there exist unique points, Xb G dQ>, yb € 90*, 
such that Xm = Vb on the c-segment, i Xb , joining x to Xb, with respect to yb, (except at the 
endpoint Xb ) with the resultant mappings being C 2 diffeomorphisms from dVt r to <90, <9fi* 
respectively, for any r < 5. 

Lemma 5.1 will follow from Lemma 7.1 by modification of u\ outside O and mollifica- 
tion. To prove Lemma 7.1, we first take any c-function in the set S, with c-normal image 
y G O* — u* and increase it until its graph meets that of u$ on O at a point Xb, which 
will lie in <9Q, since Uq is uniformly c-convex. Accordingly we obtain a c-function h G S, 
given by 

(7.5) h(x) = h Xbj y(x) = c(x, y) - c(x b , y) + u (x b ). 

Since h < u$ in O and h(xb) = uo(xb), we see that the point y must lie on the c*-segment 
£* b ,with respect to Xb, starting at yo,b = Tuo(xb) and given by 

(7.6) c x (x 6 , t Xb ) = {Du (xb) + t>y(x b ) : t > 0}. 
Conversely, for any x b G dQ : y G £* b , we have 

(7.7) h < uq in O, 

by virtue of (7.3), (taking x$ = Xb.yo = yo,b)- This proves that u\ is indeed a c-convex 
extension of u to O 5 . 

To proceed further, we let yb be the unique point in 90*, where £* b intersects 90*. 

Since u* is also uniformly c*-convex, £* b only intersects du* at the initial point yo,6- 

Actually, only the uniform c-convexity of uq is needed to justify this. Henceforth we 

restrict £* b to the segment joining y ,6 to y&. From our argument above, the mapping 
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from x b to y b is onto 90*. We claim it is also one-to-one. For suppose the c- function 
h = h XbjVb meets 90 at another point x' . By increasing h, we infer the existence of a 
further boundary point x" ', which is a saddle point for uo, contradicting (7.3). It follows 
then that the mapping from x b to y b is a C 3 diffeomorphism from 90 to 90*. Next, if 
B r is a sufficiently small tangent ball of O at £5, it will also be uniformly c-convex so 
again by (7.3), we obtain 

(7.8) h Xbiyb (x) > h Xbiy (x) VxeB r ,yet Xb . 
and thus we have 

(7.9) ui = max{no,/i Xb ,yJ, x E Q s . 

To complete the proof of Lemma 7.1, we need to show that for each x E O* 5 — O, there 
exists a unique x b E 90 where the maximum in (7.8) is attained. For this purpose, we 
invoke the c-transform of u±, 

(7.10) v (y) = sup{c(x,y) - Ul (x) : x E O* 5 }, y e 0*, 

which extends the c-transform of Wo in cu*. Moreover, by (7.7), we see that for y e 
the sup is attained at Xb- Hence 

(7.11) v (y) = c(x b ,y) - u (x b ) M y e t Xb . 

One easily verifies that vq is smooth in O* — duo* . Using (7.11) and arguing as before, 
we infer that for any point x E O" 5 — O , there exists a unique point y b E 90* such that 

K, yb (y) = c ( x , y) - c ( x , Vb) + Mvb) 

(7.12) <v (y) VyGO*. 
Moreover x lies on the c-segment, t Vb given by 

(7.13) c y (£ yb ,y b ) = {Dv (y b )+t>y*(y b ) : t E [0,5]}, 

where 7* denotes the unit outer normal to 90* and 5 is a small constant. Note that 
Dvo(y b ) = c y (x b , y b ) From (7.12), we see that the maximum in (7.9) is attained at x b , y b 
so 

(7.14) u^x) = c(x,y b ) - c(x b ,y b ) + u (x b ) x E £ yb , 

with Xuii^yt — {xb}) = yb, XuA^-b) = £% b - From the obliqueness of l Vb on 90, we have 
that the mapping from x E O r to x b is one-to-one for sufficiently small r. This completes 
the proof of Lemma 7.1. 
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From Lemma 7.1, we see that the function u\ is smooth in VI s — dfl. Furthermore, 
with 5 sufficiently small, u\ will be tangentially uniformly convex on <9fi r , that is 



(7.15) [Dijui - Cij(-,Tui)]TiTj > A , 

where r is the unit tangent vector on dfl r and Ao a positive constant. To take care of 
the normal direction, we modify u\ in O 5 — O, by setting 

(7.16) u = u 1 + bd 2 , 

where b is a positive constant and d denotes distance from O. Again for 6 sufficiently 
small, we infer that u satisfies 

(7.16) [Diju - <•;/(•. Tu)]^ > A , 

in Q s — 90 for a further positive constant Ao, and any unit vector £. 

We complete the proof of Lemma 5.1 by mollification. Let p G Cq°(Bi(0)) be a 
mollifier, namely p is a smooth, nonnegative, and radially symmetric function supported 
in the unit ball -Bi(O) such that the integral f Bl ^ P = 1- We show that a mollification 
of tt, given by 

(7.17) u e (x) = p * u 

= [ e- n p{^-)u{y)dy 

= / p(y)u(x - ey) dy 

is uniformly c-convex in O 5 / 2 , provided e < \5 is sufficiently small and x G O 5 / 2 (so that 
the value of u outside O 5 is irrelevant). Note that the image of the c-normal mapping 
of u s in is a smooth perturbation of O*, and so is also uniformly c-convex provided 
e > is sufficiently small. 

It is easy to verify that 

(7.18) Du e {x)= I p(y)Du(x-ey)dy, 

(7.19) 

D 2 u £ (x)= [ p(y)D 2 u(x-ey)dy+ [ ^p(^-) 7 ' (D+u - D~ u)(y) 
> / p(y)D 2 u(x - ey)dy, 
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where 



D + u(y) = lim Du(y'), 

D~u(y) = lim Du(y'). 

y'€il,y'->y 

Since, u* dO*, we have 

(7.20) D+u - D~u > C > on 90 

for some positive constant Co. We divide O" 5 / 2 into three parts: O" 5 / 2 = U\ U U2 U C/3, 
where 



C/i = {a; G 5/2 

U 3 = {x e n 5/2 



dist(x, 90) > e}, 
dist(;r,90) G ((l-<r)e,e)}, 
dist(o;,90) < e'}, 



where a G (|, 1) is a constant close to 1. Since u is smooth, uniformly c-convex away 
from 90, u £ is obviously smooth, uniformly c-convex in U\ provided e is sufficiently 
small. By taking a > sufficiently close to f, for any pointao G L/2, -D , u e (xo) is a small 
perturbation of Du(xq)- By (7.f9), we also see that the matrix 



(7.21) {D 2 u £ (x ) - A(x , Du e (x ))} > 0, 

namely u e is smooth, uniformly c-convex in Ui- 

Finally we verify (7.21) in U 3 . For any point x G Us, we choose a coordinate system 
such that xq = (0, • • • , 0, x , n ) 5 the origin G 90 and 90 is tangent to {x n = 0}. To 
verify (25), we first consider a tangential direction r, namely r is a unit vector tangential 
to 90 at 0. Without loss of generality we assume that r = (1, 0, • • • ,0). Then we need 
to prove that 

(7.22) D llU £ (x ) - Anixo, Du e (x )) > 0. 

By our choice of coordinates, Diu £ (xq) is a small perturbation of D\u{xq). Hence it 
suffices to verify that 

(7.23) D llU £ (x ) - An(xo, D lU (x ), D'u e {x )) > 0, 

where D'u e = (D 2 u £ ,--- ,D n u s ). By A3w, An is convex in D'u 6 , whence it follows 
readily that 

A 11 (xo,D 1 u(xo),D'u £ (x )) < / — p( )A 11 (x ,D 1 u(x ),D'u(y))dy. 
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Inequality (7.22) now follows from (7.19) and (7.23). Note that the argument also applies 
to any direction r\ provided r\ ■ 7 is sufficiently small. Next we observe from the second 
integral in (7.19) that (7.21) holds in the normal direction 7 = e n . Furthermore, 

(7.24) D nn u £ (x ) - A nn (x , Du £ (x )) > K 

for some K as large as we want, provided e is sufficiently small. Now suppose the 
least eigenvalue of the matrix (7.21) is achieved in direction £. We can decompose 
£ = C\T + C2e n . If C2 > Co for some constant cq > 0, then the matrix (7.21) in direction 
£ is positive by (7.24). Otherwise the proof of (7.22) applies and we also see that the 
matrix (7.21) in direction £ is positive. By appropriate adjustment of O, we complete 
the proof of Lemma 5.1 and consequently also Theorems 1.2 and 1.3. 

To conclude this section we show that Lemma 5.1 may be proved independently of 
the arguments in Section 5 by direct construction of a uniformly c-convex function, uq. 
To do this we let yo be a point in fi* and uq be the c*- transform of the function 

(7.25) ^y) = -(r 2 -\y-yo\ 2 ) 1/2 , 

given by 

u Q {x) = sup{c(x, y) - ip(y), y G B r (y )}, 

for sufficiently small r > 0. Then uq is a locally uniformly c-convex function defined in 
some ball Br(0), with R — > 00 as r — > 0, and the image of its c-normal mapping, 

u* :=Tu (n) C B r (y ), 

where Tuq is a diffeomorphism between O and u*. As Tuq is defined on the ball -B_r(0) D 
O, iv* is a smooth domain. Locally uq is a smooth perturbation of the c-function 

(7.26) h (-) = c(;y ) + a , 
for some constant do- 

8. Examples 

We repeat and expand somewhat the examples in [MTW], taking account that our 
cost functions are the negatives of those there. 

Example 1. 

(8.1) c(x, y) = -a/1 + \x - y\ 2 
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Here the vector field Y and matrix A are given by 



(8.2) Y(x,p) = x+ !> 



1 — |p 



2' 



A(x, p) = A(p) = - (1 - \p\ 2 ) 1/2 (/ - p <g> p) . 

The cost function satisfies condition A3. We remark that condition Al is only satisfied 
for \p\ < 1 but this does not prohibit application of our results as the boundedness of 
target domain fi* ensures that \Du\ < 1 for solutions of (1.7), (1.14). More generally the 
conditions p,q G M, n in Al may be replaced by p, q G some convex sub-domain. 

Example 2. 

(8.3) c(x,y) = -^l-\x-y\ 2 

Here c is only defined for \x — y\ < 1. The vector field Y and matrix A are given by 

(8.4) Y(x,p) = -x+ 7=== , 

V 1 + \P\ 

A(x, p) = A(p) = (l + \p\ 2 ) 1/2 (J + p ® p) . 

The cost function satisfies condition A3. In order to directly apply our results we need 
to assume O and are strictly contained in a ball of radius 1. 



Example 3. Let /, g G C 2 (fi),C 2 (fi*) respectively and 
(8.5) c(x, y) =x-y + f(x)g(y) . 

If |V/.V<7| < 1, then c satisfies Al, A2. If /, g are convex, then c satisfies A3w, while if 
f,g are uniformly convex, then c satisfies A3. As indicated in [MTW], the function (8.5) 
is equivalent to the square of the distance between points on the graphs of / and g. 



Example 4. Power costs. 

(8.6) c(x,y) = ± — \x — y\ m , m^O, log\x — y\ , m = 0). 

For m/1 and x ^ y, when m < 1, the vector fields Y and matrices A are given by 

(8.7) y(x,p) = x ± |p| m - 1 p, 

A(x,p) = A(p) = ±{ |p| ^-i / + (m — 2) \p\~ m-ip ®p|. 
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The only cases for which condition A3w is satisfied are m = 2(±) and — 2 < m < 1 (+ 
only ). For the latter, condition A3 holds for — 2 < m < 1. To apply our results directly 
in the latter cases, we need to assume O and O* are disjoint. 

In [MTW] we also considered the cost function 

(8.8) c(x,y) = -(1 + \x-y\ 2 ) p/2 

for 1 < p < 2, extending Example 1 to p > 1. We point out here that these functions only 
satisfy A3 under the restriction \x — y\ 2 < ^-j-. This condition was omitted in [MTW]. 

Example 5. Reflector antenna problem 

Corresponding results and examples may be obtained on other manifolds such as the 
spheres S n . Indeed the considerations in [MTW] stemmed from the treatment of the 
reflector antenna problem by Wang in [Wl], which may be represented as an optimal 
transportation problem on the sphere S n with cost function 

(8.9) c(x,y) = log(l-x-y), 

which is simply the spherical analogue of the case m = in Example 4 above. The 
corresponding vector field Y is now given by 

(8.10) Y(x,p) =x- T-^ix + p), 

1 + \p\ z 

where now p belongs to the tangent space of S n at x, while the matrix A is given by 

(8.11) A = ^(\p\ 2 -l)g -pxp, 

where go denotes the metric on S n . See [Wl, W2,GW, MTW] for more details. When the 
domains Q and Q* have disjoint closures, and spherically uniformly convex boundaries, 
we obtain the global regularity of potentials. 

We will defer further examination and extensions to intersecting domains and other 
cost functions in a future work. We also point out here that Example 4 provides regularity 
for quadratic cost functions on spheres when the points x and y are sufficiently close. 
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